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Abstract 



An investigation of CP violation was performed using a total of 24 candidates for 
B° — > J/^Kg decay, with a purity of about 60%. These events were selected from 4.4 
million hadronic Z° decays recorded by the OPAL detector at LEP. An analysis proce- 
dure, involving techniques to reconstruct the proper decay times and tag the produced 
b-flavours, B° or B°, has been developed to allow a first direct study of the time dependent 
CP asymmetry that, in the Standard Model, is sin 2(3. The result is 

sin 2/3 = 3.2^ ±0.5 , 

where the first error is statistical and the second systematic. This result is used to 
determine probabilities for different values of sin 2(3 in the physical region from —1 to +1. 
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1 Introduction 



CP violation was observed in 1964 in K° decays |IJ, and, so far, this phenomenon has been seen 
only in the K system. CP violation can be accommodated in the Standard Model, provided 
that the CKM matrix elements are allowed to be complex. The CP violating effects associated 
with b hadrons are expected to be larger than in the K system 0. It is therefore important 
to investigate CP violation in the B system, where the relation between CKM matrix elements 
and CP violation can be tested. Previous studies, yielding null results || £§], have focussed on 
CP violation in inclusive B decays, predicted to be at the level of 10 -3 . 

This paper presents a study of CP violation in B° — > J/^Kg decays. The decay mode J/ipK^ 
has long been considered a 'golden' channel for CP violation studies [|[ |J, since the final state 
is a CP eigenstate which is experimentally favourable for reconstruction because the J/ip and 
Kg are narrow states and J/ip —>■ £ + £~ decays give a distinctive signature. In addition, CP 
violation in this channel is dominated by diagrams having a single relative phase, allowing a 
clean extraction of the phase of a CKM matrix element. In the Standard Model, the expected 
time-dependent rate asymmetry, A(t), is given by 

B°(t) - B°(t) 

A{t) = B o(t) + BO(t) = - sin2 ^ sinA ^ , (1) 

where the parameter Awid is the mass difference between the two B° mass eigenstates and B°(£) 
(B°(t)) represents the rate of produced B°'s ( B°'s ) decaying to J/ipK^ at a given proper decay 
time, t. The angle ft is given by 



ft = arg 



-argVt d . (2) 



Constraints on the CKM matrix, including measurements of CP violation in the K system, 
imply that the Standard Model expectation for sin 2ft lies in the range 0.3-0.9 0. Other 
models of CP violation, such as the Superweak model @, would give a time dependence of the 
same form, but with sin 2ft replaced by another amplitude of magnitude less than or equal to 
one. 

At LEP, due to the small branching ratios of B° — ► J/ipK^ and J/ip — > £ + £~ , only a handful 
of these decays may be seen. It is therefore important to maximise the reconstruction efficiency 
and to determine the b-flavour at production with a minimum error rate. In contrast, the proper 
time resolution is not critical, since the frequency of B° oscillation is easily resolved. In this 
analysis, B° — > J/^Kg decays are reconstructed and their decay proper times are measured. 
The production flavour (B° or B°) of each candidate is determined using a combination of 
jet and vertex charge techniques. The CP- violating amplitude, sin 2ft, is extracted using an 
unbinned maximum likelihood fit to the proper-time distribution of the selected data. Such a 
fit has higher expected sensitivity than a time-integrated fit, even if the domain of integration 
of the latter is optimised, since it uses the additional, well measured, decay-time information. 

The next section describes the event selection and proper time estimation. Section 3 de- 
scribes the tagging of the production flavour. In section 4, the fits and results are presented, 
with systematic errors discussed in section 5. Discussion of the result and conclusions are given 
in section 6. 
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2 B° — > J/^Kg reconstruction 



A detailed description of the OPAL detector may be found elsewhere f7j. The basic selection 
of B° — > J/^Kg decays was described in a previous publication ||, which used data collected 
between 1990 and 1994 to identify various exclusive B decay modes. For this analysis, the 
1995 data were included to give a total of 4.4 million events passing the basic hadronic event 
selection. The efficiency of the B° — ► J/ipKg selection was also improved (and applied to the 
full data sample) by relaxing or modifying the criteria, at the expense of a somewhat larger 
background. The efficiency and purity of the selection was studied using Monte Carlo events 



generated with Jetset 7.4 |§ and processed through the OPAL detector simulation |l(J 



Lepton candidates were required to satisfy the polar angleQ cut | cos#| < 0.97 and to have 
track momenta larger than 1.5 GeV/c (2 GeV/c ) for muon (electron) candidates. Muons were 
identified by requiring an extrapolated track to match the position of a muon chamber segment 
to within 4 standard deviations. Muon candidates were also considered if the muon chamber 
segment had no ^-coordinate reconstructed, provided that the matching in position and angle 
were within 4 standard deviations in the r-0 plane. As in ||, muons identified in the hadronic 
calorimeter were accepted in regions without muon chamber coverage. 



Electrons were identified [D] and photon conversions rejected |12| using artificial neural 
network algorithms. When the electron energies are determined only from the reconstructed 
track momenta, photon radiation causes the reconstructed J/ip mass spectrum to have a long 
tail to lower masses, reducing the efficiency of a mass cut. Therefore the electron energies 
were determined using in addition the information from the lead-glass calorimeter. The energy 
contained in a cone of 30 mrad around the impact point of the electron track on the calorimeter, 
plus energy contained in lead-glass blocks touching this cone, were summed Jl3| . This sum was 



used as the energy of the electron if larger than the track momentum, otherwise the track 
momentum was used. 



J /ijj candidates were selected by demanding two electron or two muon candidates of opposite 
charge with an opening angle of less than 90°. The invariant mass of the two leptons was 
required to lie in the range 2.95-3.25 GeV/c 2 for J/ip — > fi + pT candidates and 2.95-3.40 GeV/c 2 
for J /if) — > e + e~ candidates (in the latter case over-correction for photon radiation causes a 
significant tail at higher masses). 

The K s selection was based on the procedure described previously jnj, considering the 
intersection of all track pairs of opposite charge (excepting J/ip candidate tracks) passing certain 
quality criteria. The projection of the the K s momentum vector in the r-0 plane was required to 
point back to the beam spot position to within 8°. The beam spot position was measured using 
charged tracks from many consecutive events, thus following any significant shifts in beam 
position during a LEP fill fll6| . The impact parameter^ significance (the impact parameter 
divided by its error) in the r-0 plane of the Kg with respect to the J/tp vertex was required to 
be less than 5. The reconstructed distance between the J/ip vertex and the Kg decay vertex, 

1 The right-handed coordinate system is defined with positive z along the e _ beam direction, x pointing to 
the centre of the LEP ring, 6 and <f> as the polar and azimuthal angles, and r 2 = x 2 + y 2 . The origin is taken 
to be the centre of the detector. 

2 The impact parameter of a track with respect to a vertex is defined as the distance of closest approach of 
the track to that vertex. 
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divided by its error, was required to exceed 2. If the Kg decay vertex was inside the active 
volume of the jet chamber (r > 30 cm), the radial coordinate of the first jet chamber hit on 
either track was required to be at most 10 cm less than the decay vertex radial coordinate, 
and the tracks were required not to have any associated vertex chamber or silicon microvertex 
detector hits. The invariant mass of the Kg candidate was required to lie in the range 0.45- 
0.55GeV/c 2 . In order to suppress a potential background from Ab — > J/if)A decays, the Kg 
candidate was rejected if its invariant mass under either proton-pion hypothesis was in the 
range 1.110-1.121 GeV/c 2 . 

B° candidates were reconstructed by combining J /if) and Kg candidates from the same 
thrust hemisphere^. Kinematic fitting using the SQUAW package |TJ]], constraining the J /if) 
and K s masses to their nominal values, was employed, and the probability of the kinematic 
fit was required to exceed 1%. The invariant masses of selected B° candidates were required 
to lie in the range 5.0-5.6 GeV/c 2 , and the B° energies to exceed 20GeV. The efficiency for 
reconstructing the decay B° — ► J/^K S — > £ + £~tt + tt~ was estimated to be 17.3 ± 1.4% where the 
error is due to Monte Carlo statistics. This compares to 10.9% for the previous selection ||. 

The distribution of reconstructed mass is shown in Figure |l](a) for candidates passing the 
entire selection except the J/ipK^ mass cut. In total, 24 candidates were selected in the mass 
region 5.0-5.6 GeV/c 2 . The background was estimated from the data using an unbinned max- 
imum likelihood fit to the joint distribution of the J/^Kg mass and energy, including mass 
values between 4 and 7GeV/c 2 . The shape of the mass distribution for the signal was taken 
from Monte Carlo, and parametrised using two Gaussians for the peak of the distribution and 
a third to account for a significant tail to lower masses. The shape of the background mass 
distribution was also taken from Monte Carlo, and parametrised using a polynomial function. 

The mass and energy distributions were taken to be uncorrelated both for the signal and 
the background, since the correlations seen in Monte Carlo were small and had a negligible 



effect on the fit. Peterson fragmentation functions |L7| were used to parametrise the energy 
distributions. For the signal, the Peterson parameter e was taken from a fit to the Monte Carlo 
energy distribution, while for the background it was allowed to vary. In total, three parameters 
were allowed to vary in the fit: the number of signal candidates, the position of the B° mass 
peak and the Peterson parameter for the background energy distribution. The result of the fit 
is shown in Figure |l|(a) as a function of mass, in Figure |l](b) as a function of energy for masses 
between 5.0 and 5.6 GeV/c 2 (signal region) and in Figure |TJ(c) as a function of energy for masses 
outside this region (sideband). The overall fitted purity of the 24 B° — > J/if>Kg candidates is 
60 ±8%, where the error is statistical. The fit was used not only to determine the overall purity 
of the sample, but also to assign event-by-event background probabilities, /bac, to be used in 
the fit for sin 2/3, according to the reconstructed J/^Kg mass and energy of each candidate. 

The B° decay length, /bo, was determined from the distance between the average beam 
spot position and the J /if) vertex in the x-y projection, correcting for the polar angle using 
the direction of the J/ipK^ momentum vector. The B° momentum, p#o, was taken from the 
constrained fit of the J/^Kg system, and the proper decay time was then calculated as 

tree — Ib° ' j (3) 

Pbo 

The two hemispheres were separated by the plane perpendicular to the thrust axis of the event and 
containing the e + e _ interaction point. 



6 



i — i — i — i — I — i 1 — i 1 — I — i — i — i 1 — I 1 — i — i — i — I — i — i — i r 




+ OPAL data 

o o 
m B%JAl/K s 

^ B%JA|/K*° 

□ background 



6 6.5 7 

2 

Invariant mass (GeV/c ) 



c — I — I — I — I — I — I — I — I — I — I — I — I — I — I — q 



Signal 
region 




> 

o 

IT) 



a 



40 50 

Energy (GeV) 



35 
30 



25 



W 20 



15 
10 
5 




~i — i — i — i — i — i — i — i — i — i — i — i — i — r 



(C) 



Sideband 
region 



i i i i 



20 



30 40 50 

Energy (GeV) 



Figure 1: (a) The mass distribution of J/^Kg candidates, (b) The energy distribution of 
the candidates with masses between 5.0 and 5.6GeV/c 2 . (c) The energy distribution of the 
candidates with masses outside this region. In each case, the data are shown by the points 
with error bars, and the fit is shown by the open histogram. The estimated contribution from 
the B° — > J/V>K S signal is shown by the cross hatched histogram. The estimated contribution 
from B° — > J/ipK*° is also shown in (a) by the hatched histogram. Note that the J/ipK*° 
contribution is taken directly from Monte Carlo, while the shapes of the other distributions are 
parametrised (see text). 
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where M B o is the B° mass, taken to be 5.279 GeV/c [ 18| . The uncertainty on t rec , typically 



o t — 0.1 ps, was estimated by combining the uncertainties on Z B o and p-^o. 



3 Tagging the produced b-flavour 

Information from the rest of the event was used to determine the production flavour of each 
candidate. The weighted track charge sum, or 'jet charge', gives information on the charge, and 
hence b-flavour, of the primary quark initiating the jet within which the J/ipK^ candidate was 
isolated. Additionally, since the Z° decays into quark-antiquark pairs, measuring the b-flavour 
of the other b quark produced in the event also provides information on the production flavour. 

In this analysis, three different pieces of information were used to determine the B° pro- 
duction flavour^: (a) the jet charge of the highest energy jet other than that containing the 
B° candidate, (b) the jet charge of the jet containing the B° candidate, excluding the tracks 
from the J/if> and K s decays, and (c) the vertex charge of a significantly separated vertex (if 
existing) in the opposite hemisphere. 

Jets were reconstructed from tracks and electromagnetic clusters not associated to tracks 



using the JADE E0 recombination scheme [[19] with a y cu t value of 0.04 



The jet charge for the highest energy jet other than the one containing the J/^Kg candidate 
was calculated as : 

Qopp = E(VM *, (4) 

j \ -'-'beam / 

where p\ is the longitudinal component of the momentum of track i with respect to the jet 
axis, -Ebeam is the beam energy, g« is the electric charge (±1) of each track and the sum is 
taken over all tracks in the jet. Using simulated data, the value of k that minimised the mistag 
probability (the probability to tag a B° as B° and vice versa), was found to be approximately 
0.5. The mistag probability includes the effect of B° mixing in this jet. The jet charge for 
the jet containing the B° candidate, Q sam e, wa s calculated in the same way, but excluding the 
J /ifj and Kg decay products. These particles contain no information on whether their parent 
was produced as a B° or B°, and would only dilute the information from the fragmentation 
tracks. The optimal value of k was found to be 0.4 in this case, smaller than that for the other 
b-hadron as the high momentum B decay products were excluded. 

Secondary vertices were reconstructed in jets in the hemisphere opposite to the B° candidate 
in data where silicon microvertex information was available, using the algorithm described in 



20[ . For the 1991 and 1992 data, vertices were reconstructed in the x-y plane only. In the 
1993-1995 data, precise z coordinate information from the silicon microvertex detector was also 
available, and vertices were reconstructed in three dimensions using an extension of the vertex 



finding algorithm as in |2T| . A secondary vertex was accepted if the distance from the primary 
to the secondary vertex divided by its error (the vertex significance) exceeded 3. If more 
than one vertex in the opposite hemisphere satisfied this requirement the one with the highest 



4 Leptons in the opposite hemisphere were also investigated as a possible tag, but in this sample the events 
with selected leptons were found to have large jet charges (in agreement with the lepton charge), so that there 
was no significant gain in tag purity from the use of leptons. 
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significance was taken. Approximately 40% of Monte Carlo B° — > J/^Kg events had such an 
accepted secondary vertex in the opposite hemisphere. For the selected vertex, the charge Q vtx 
and its uncertainty <jQ vtx were calculated using an improved version of the algorithm described 
in p2 |. For each track i in the jet containing the vertex, the track momentum, the momentum 
transverse to the jet axis, and the track impact parameters with respect to the primary and 
secondary vertices in the r-(f) and r-z planes (the r-z information was only used in 1993-1995 
data) were combined to form a weight Wi using an artificial neural network algorithm. The 
weight Wi quantifies the probability for track i to belong to the selected secondary vertex. The 
vertex charge is then calculated as: 

Qvtx = W iQi ' ( 5 ) 
i 

and the uncertainty as: 

^Qvtx = -Wi)q 2 i • (6) 



For events with such a selected secondary vertex (9 of the 24 B° candidates fall into this 
category), a neural network was constructed to tag the produced B° or B°, combining the four 
inputs, Q v tx, o"Qvtx, Qojip ' 5 an d Qtlme- The network was trained on a large Monte Carlo sample, 
and a variable 

Qb(x) = f B o(x) - /go (x) (7) 

was calculated as a function of the network output, x, where /b ^) (fB°( x )) * s the fraction 
of candidates at a particular value of x due to produced B° (B°) according to Monte Carlo 
(which included the effects of B° mixing). The variable Qb represents the effective produced 
b-flavour for each candidate (Qb = +1 or —1 for pure B° or B°, respectively), and \Qb\ — 1 — 2r? 
the tagging dilution, where rj is the probability to tag the production flavour incorrectly. The 
average value of \Qb\ is 0.38 for such events, according to Monte Carlo. 

For events without such a selected secondary vertex, only the jet charge information is 
available. In this case the two jet charges were combined linearly to form 

Q2jet = Qsaine ~ 1-43 • Qop P °' 5 ) (8) 

where the factor of 1.43 was optimised using simulated data to minimise rj. Distributions 
/ B °(<52jet) and /B°( ( 52jct) were formed and QB(Q2jct) was determined by analogy to equation [7]. 
For these events, the Monte Carlo predicts the average value of \Q B \ to be 0.31. 

It is important to ensure that the tagging dilution arising from the jet and vertex charges, 
which are used to determine the Qb values, is correctly modelled by the Monte Carlo. The 
tagging dilutions arising from Qopjp ' 5 an d from Q v tx were checked using large samples of data 
and Monte Carlo inclusive lepton events, selected as in [|J. In such events, the charge of 
the lepton (usually from a semileptonic decay of a b hadron), is strongly correlated with the 
produced b-flavour. The distributions of <5opp°' 5 an d Qvtx, multiplied by the lepton charge Qt, 
were compared for data and Monte Carlo, and found to be consistent (see Figures 0(a) and (b)). 
The mean values of these distributions (which are sensitive to the tagging dilution of Q^ ' 5 or 
Qvtx) for data divided by those in the Monte Carlo, were found to be 1.13 ±0.09 and 1.05 ±0.11 
for Qq™ ' 5 and Q vtx respectively. The Q^Zme dilution was checked using samples of B° — > D* ± £ 



candidates selected in data and Monte Carlo, as in p3]. (In calculating Q^e 4 , tracks from 
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the D* ± £ combinations were excluded.) The average dilution seen in the data divided by that 
observed in Monte Carlo was found to be 0.98 ± 0.20. The distributions of Qg ■ Qsame 4 f° r data 
and Monte Carlo are shown in Figure 0(c). Only events with reconstructed B° decay time 
tree < 2 ps are included, to reduce the tagging dilution due to B° mixing. 

The samples of D* ± £ events were also used to check the tagging dilution of the Qb values 
(formed from the combinations of jet and vertex charges). The average dilution seen in the 
data divided by that seen in Monte Carlo was found to be 0.96 ± 0.14. The distributions of 
Qt ■ Qb are shown in Figure 0(d), again for events with t rec < 2 ps. Both Figures 0(c) and 0(d) 
show agreement between data and Monte Carlo. 

The Q vtx , Qopp°' 5 and Qsame 4 distributions are not necessarily charge symmetric because of 
detector effects causing a difference in the rate and reconstruction of positive and negative 
tracks. These effects are caused by the material in the detector and the Lorentz angle in the 
jet chamber. They were removed by subtracting offsets from the Q vtx , Qopp ' 5 and Q^ame 4 values 
before the Qb tagging dilutions were calculated. The Q vtx and Qopp°' 5 offsets were determined 
using the inclusive lepton events selected from data. The Q^^ 4 offset was determined from 
Monte Carlo B° jets, since no pure sample of fully reconstructed B° decays is available from 
the data. However, the D*^£ events do allow this offset to be checked, even though some 
extra tracks may be present from D** decays. The data and Monte Carlo were found to be 
in good agreement. The normalised offsets (the offsets divided by the r.m.s. widths) of the 
charge distributions were found to be +0.029 ± 0.011, +0.018 ± 0.007 and +0.036 ± 0.018 for 
Qvtx, Qopp ' 5 and Qs^ie 4 respectively The error quoted for the Q^ e 4 offset is the statistical 
precision of the D £ events in data. If these offsets were not removed, they would induce 
respective shifts of —0.004, —0.003 and +0.009 in the overall Qb distribution. If no corrections 
were applied for the offsets, the combined shift would be +0.002. 



4 Fit and results 

The reconstructed proper time t rec and tagging variable Qb of each of the 24 candidates is 
shown in Figure 0(a). The events with invariant mass in the range 5.15-5.40 GeV/c 2 are shown 
as filled circles, whilst the events in the ranges 5.00-5.15 and 5.40-5.60 GeV/c 2 , which have 
lower signal purity, are shown as open circles. 

In order to quantify the CP asymmetry in the data, an unbinned maximum likelihood fit 
was constructed, using four inputs for each B° — ► J/VKg candidate: t rec , a t , Qb and the event- 
by-event background probability /bac, derived from the mass and energy of each candidate (see 
Section 2). The total likelihood for an event was given by 

(l ,/bac) ' ^sig + /bac ' ^bac • (9) 

The signal likelihood C slg was defined as: 

£ sig (t r ec, o- t ,Q B ;sin2/3, Am d ,r°) = .F sig (t) ® G(t - t rec , a t ) , (10) 

where t is the true proper decay time, and G(t — t Tec , a t ) is a Gaussian representing the proper 
decay time resolution. The true proper time distribution is given by 

J zr sig (t;sin2/3, Am d ,r°) = ex P(~*/ T } . (1 - g B s i n 2/3sin Am d t). (11) 
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Figure 2: The distribution of the product of the lepton charge Qe and (a) the Qopp' 5 jet 
charge, (b) the Q v t x vertex charge, (c) the Q^Zme charge and (d) the Qb tagging variable. 
Distributions (a) and (b) are taken from an inclusive lepton sample, whilst distributions (c) 
and (d) are taken from a B° — > D*^£ sample, only including events with reconstructed proper 
time t rec < 2 ps. The data are shown as points and the Monte Carlo predictions as histograms. 
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The B° lifetime, r°, was taken as 1.56±0.06 ps |H|, and Am d was taken as 0.467±0.022lg;g^ ps 



-i 



[24|. The likelihood for the background, £bac> is defined in the same way, with the true proper 



time distribution: 

T v exp(-t/r bac ) 

-/"bac(* STbacJ = • (12) 

7~bac 

The background is dominated by bb events, and is assumed to have no CP asymmetry. Possible 
bias due to this assumption is treated as a systematic error. The effective background lifetime, 
Tbac, was taken to be 2.0 ps from the Monte Carlo background sample. This value is larger 



than the average b lifetime of 1.55 ps |L8j because the energy of the b hadron is systematically 
underestimated for the background events, since the tracks assigned to the J/^Kg candidate do 
not, in general, include all the b-hadron decay products, and include fragmentation products. 
A large variation of ±0.4 ps in this parameter is considered in the systematic errors. 

Fitting the data for the single parameter sin 2(3, gave the result 

sin 2/5 = 3.2tH . 

The corresponding Alog£ distribution is shown in Figure 0(b). It can be parametrised as 

-Alog£ = 0.116(sin 2/5- 3. 2) 2 + 0.00224(sin 2(3- 3.2) 4 sin2/5<3.2 
-Alog£ = 0.125(sin 2/5- 3. 2) 2 + 0.00985(sin 2(3- 3.2) 4 sin2/5>3.2. 

The parametrisation can be used to combine this result with future results from other exper- 
iments. To compare the fitted result with the data, an estimator, A, of the B° — > J/^Kg 
asymmetry (corrected for the average dilution in each time bin) is shown in Figure |3|(c) with 
the fit result superimposed, where 

A _ (1 — /bac) • Qb n „\ 
E (1 - /bac) 2 " Ql ' (13) 

and the summations are over all the events in a given time bin. The large observed values of A, 
typically exceeding the physical range of the asymmetry, are due to the tagging dilution factors 
and, to a lesser extent, the background fraction. 



5 Systematic errors and cross checks 

The main sources of systematic error and their effect on the measurement of sin 2(3 are listed in 
Table [l]. The fit result is sensitive to the level and possible CP asymmetry of the background, 
the accuracy of the decay time reconstruction and the production flavour tagging dilution. 

• The event-by-event purities of the 24 candidates have significant statistical errors from 
the background fit described in Section 2. The three fitted parameters (the number of 
signal candidates, the position of the B° mass peak and the Peterson parameter for the 
background energy distribution) were each varied by their statistical errors, one at a time, 
and the effect on sin 2(3 determined. As the correlations between these parameters were 
found to be less than 20%, these effects were added in quadrature, leading to a total error 
of on sin 2(3. 
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Figure 3: (a) The distribution of Q B versus t rcc for the J/^Kg candidates with invariant masses 
in the range 5.15-5.40 GeV/c 2 (filled circles), and 5.00-5.15 GeV/c 2 and 5.40-5.60 GeV/c 2 (open 
circles); (b) the — Alog£ value as a function of sin 2(3 from the fit to the 24 data candidates, 
with the physical region indicated by the dotted lines; (c) the distribution of the corrected 
asymmetry, A, versus t rcc for the J/?/>K s data, with the fit result superimposed. 
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The background mass distribution is taken from a fit to a Monte Carlo sample with 
four times the statistics of the data sample. Various different parametrisations and the 
binned Monte Carlo distribution itself were tried. As an alternative method, the shape 
of the background was taken to be a falling exponential, with both the normalisation and 
decay constant of the background being fitted to the data, as in ||. This fit predicts a 
signal purity of 52% compared to 60%. The Monte Carlo predicts a significant departure 
from the exponential shape due to decays of the type B — > J/^K^X, followed by K* — > 
Kg7r when K* are produced. Including this contribution explicitly and letting the fitted 
exponential describe the remaining background gives a purity of 53%. The data show a 
deficit of events in the region 4.7-4.9 GeV/c 2 , possibly indicating that the Monte Carlo 
overestimates the background in this region. A further fit to the data alone was therefore 
performed, using only the data above 4.9GeV/c 2 , with the background described by a 
falling exponential, resulting in a signal purity of 73%. The largest variation in sin 2/5 
resulting from these different parametrisations was found to be ^o^. An uncertainty of 
±0.32 was taken for the systematic error due to the background parametrisation. 

Uncertainty on the assumed mass and energy distributions for the signal also affects the 
result of the background fit. Monte Carlo events with tracking resolution degraded by 
10% were used to parametrise the signal mass distribution, resulting in a shift in the 
fitted value of sin 2(3 of 0.07. The functional form used to fit the signal mass distribution 
was changed from three Gaussian functions to two, one for the peak of the distribution 
and another for the tail. This caused a shift of 0.11 on sin 2(3. The uncertainty on the B° 
energy distribution was assessed by varying the Peterson parameter e to cause a change 
in the mean scaled energy of 0.02, larger than the uncertainty on the mean scaled energy 



of B hadrons [25|. The effect on sin 2(3 was negligible. 



The final state from the decay B° — > J/ipK*° followed by K*° — > KgV is expected to 
be mainly CP even (i.e. opposite CP to the J/V'Kg final state), and so could give rise 
to a possible CP asymmetry in the background. The contribution from such decays in 
the signal region was estimated from Monte Carlo to be 0.7 events, and is indicated 
in Figure |l|. If such a contribution had a maximal asymmetry, the effect on the fitted 
sin 2(3 would be 0.03. The contribution to the background CP asymmetry from B° decays 
involving K° mesons was found to be negligible. 

The event-by-event proper time resolution at is used in the likelihood fit. Monte Carlo 
studies indicate that the distribution of errors in reconstructed proper time divided by at 
is well described by a Gaussian with zero mean and width 1.15 ± 0.15. If the proper time 
resolution is scaled by 1.3, the resulting change in sin 2(3 is 0.01. 

The description of Qopp ' 5 an d Qvtx by the Monte Carlo was tested by comparing the 
correlation of these charges with the lepton charge in inclusive lepton events in data and 
Monte Carlo, as described in Section 3. The uncertainty on the Qb values was assessed 
by scaling the Q^pp ' 5 an d Q v tx values independently by 1.16 and 1.12, respectively. These 
scalings correspond to the sum in quadrature of the differences seen between data and 
Monte Carlo and the statistical precision of the comparisons. The uncertainty on the 
modelling of Qs^ e 4 was assessed using D^i data as described in Section 3. The systematic 
uncertainty was determined by scaling the values of Q^ e 4 by 1.2, again corresponding 
to the quadrature sum of the difference seen between data and Monte Carlo and the 
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Source 


6(sin2/3) 


Background level (data statistics) 


+0.06 
-0.07 


Background shape 


±0.32 


Signal shape 


±0.13 


Background asymmetry 


±0.03 


Proper time reconstruction 


±0.01 


Jet and vertex charge modelling 


+0.31 
-0.26 


Jet charge offsets 


+0.14 
-0.08 


Vertex charge performance 


±0.01 


Amd value 


±0.10 


B° lifetime 


±0.01 


Background lifetime (±0.4 ps) 


+0.01 
-0.02 


Total 


+0.50 
-0.46 



Table 1: Sources of systematic error in the measurement of sin 2/3. 

statistical precision of the comparison. The total systematic error on sin 2(3 from these 
effects is tlH- 

• The offsets applied to Q vtx and Qopp ' 5 were determined from data as described in Section 3. 
These were varied by their statistical uncertainties. The offset to the Q^me^ J et charge 
was determined from Monte Carlo, and checked using the B° -> D*±£ candidates. The 
offset was varied by the statistical precision of this test. The effect of these variations 
results in changes in sin 2(3 of Taos- 

• The performance of the vertex charge algorithm is sensitive to the tracking resolution. 
The Monte Carlo has been tuned to reproduce the data impact parameter resolutions as 
a function of cos#, p and the different sub-detectors contributing to a track measurement. 
Residual uncertainties were estimated by degrading the resolution of all tracks by 10% 
using a simple smearing technique. The neural network training and mistag parametri- 
sations were repeated on this degraded sample, which was then used to derive the Qb 
values that enter the fit for sin 2(3 . The resulting change in sin 2(3 was 0.01. 

• The values for Am d and the B° lifetime were varied within their errors to give the un- 
certainties listed in the table. The value of Tb ac was varied by a conservative 0.4 ps (the 
difference between the predicted Monte Carlo background lifetime and the average B 
meson lifetime) to allow for uncertainties on the B energy mismeasurement. 

The total systematic error is thus ±0.5. Many of the sources of error have a statistical 
component, and many of them scale with the fitted value of sin 2(3. The systematic error would 
thus decrease in an analysis with higher statistics. 

A number of consistency checks were also performed. The result was found to be stable 
when the least well tagged events (those with \Qb\ < 0.25), the events with highest background 
(/bac > 0.5), or the events outside the purest mass region (5.15-5.40 GeV/c 2 ) were removed. 
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The values of sin 2(3 resulting from these checks were found to be 4.0±2i3 j 3. 21^0 an d 2.812TO 
respectively. The data were fitted for the B° lifetime, giving a result of 1.2±o;| ps (independent 
of sin 2/5), consistent with the world average. In addition, the assumption that the background 
exhibits no CP asymmetry was tested by repeating the fit, using only events in the sideband 
region and setting /bac to for every event. The fitted value of sin 2(3 in this case was 0. 381^49, 
consistent with zero. The selection cuts were loosened to give a sideband data sample three 
times larger, and a fitted value of sin 2(3 of —0.12 ± 0.30 was obtained. The background and 
CP asymmetry fits were also repeated on a Monte Carlo sample with no CP violation and four 
times the data statistics, giving asymmetries consistent with zero for both signal and sideband 
regions. The B° — > D £ sample was used to perform a further cross check for the absence of 
large biases in the determination of Qb- The events were separated according to the sign of Qi, 
and the average Qb, (Qb) an d (Qb), was calculated for each subsample. To account for charge 
biases due to the D* 1 * 1 selection, the average value of Qb was calculated as ((Qb) + (Qb))/2 
and was found to be —0.016 ± 0.011, consistent with zero. The Monte Carlo prediction for the 
Qopp' 5 offset, which is not used in the analysis, disagrees with the value determined from the 
data. If the fit is repeated taking all offsets from the Monte Carlo, sin 2(3 is shifted by —0.28. 
This discrepancy between data and Monte Carlo does not affect the description of the tagging 
dilution. 

The value of sin 2(3 can also be estimated from the time-integrated asymmetry. In this 
case, the lower limit of the time integration can be varied to optimise the sensitivity — i.e. the 
ability to distinguish different true values of sin 2(3. For data samples of this size and purity, 
the optimum lower boundQ was found, using Monte Carlo studies, to be 0.7 ps. The value of 
sin 2(3 obtained from our data sample using this method is 2.0^^ 5, where 12 events are included 
in the range of integration. The probability of obtaining time dependent and time integrated 
measurements disagreeing at this level or more was found to be 20%. Monte Carlo studies 
indicate that the errors obtained from both types of fit increase as the central values deviate 
from zero, and that the time dependent fit yields smaller errors on average. They also show 
that the time dependent fit has a greater sensitivity to the true value of sin 2(3 than the time 
integrated method, even after optimising the lower time-integration bound. 



6 Discussion and conclusion 



The result from this analysis can be interpreted by calculating the probabilities to see a devi- 
ation, in the positive sin 2(3 direction, at least as far from the true value as that observed, for 
different true values of sin 2(3. The deviation is defined by the difference in log£ between the 
fitted value and the assumed true value. This definition is used because the sensitivity varies 
from experiment to experiment. Monte Carlo samples of 24 candidates, with the same back- 
ground and tagging distributions as those expected in the data, were generated to determine 
these probabilities. The probabilities for the log £ differences seen in the data, with correction 
for the systematic error, were found to be be 1.6%, 7.8% and 21.3% if the true value of sin 2(3 
were —1,0 and +1, respectively. These probabilities indicate the consistency of the result with 
these values of sin 2(3, and should not be interpreted as confidence levels. The distributions 
of fitted sin 2(3 for Monte Carlo experiments with true sin 2(3 of —1, and +1 are shown in 

5 The value is smaller than that which would be obtained in the absence of background. 
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Figures |](a), (b) and (c). 



An alternative interpretation is given by the Bayesian approach assuming equal a priori 
probabilities for every possible true value of sin 2(3. In this case, the probabilities for sin 2/5 to 
be greater or less than zero, with correction for the systematic error, are found to be 68.5% 
and 31.5% respectively. 

In conclusion, the time dependent CP asymmetry in the decays B° — > J/^Kg and B° — > 
J/i/jK^ has been measured using data collected with the OPAL detector at LEP between 1990 
and 1995. From 24 reconstructed B° — > J/ipK® candidates with a purity of about 60%, the CP 
violation amplitude, which is sin 2/3 in the Standard Model, has been found to be: 

sin 2(3 = 3.2±^ ±0.5 , 

where the first error is statistical and the second systematic. The systematic error has a large 
statistical component, and much of it scales with the central value. 

This is the first direct study of the CP asymmetry in the B° — > J/^Kg system. It can be 
combined with other results in the future by using the log-likelihood curve given in section 4. 
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Figure 4: The distributions of fitted sin 2(3 in Monte Carlo samples of 24 events. The dotted 
lines indicate the data fitted value. 
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